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<^ : Abstract 



^ ■ We calculate the ground-state energy for the nucleus ^He with V4 nu- 

(•^ . cleon interactions, making use of a Jastrow description of the corresponding 

, wavefunction with state- dependent correlation factors. The effect related to 

QQ I the state dependence of the correlation is quite important, lowering the upper 

0^ ' bound for the ground-state energy by some 2 MeV. 



The physical problem of describing many interacting identical particles from a micro- 
scopic point of view can be attacked using a number of techniques. A frequent ansatz for 
the ground state (g.s.) of a many-body system is the Jastrow wavefunction. For a nuclear 
5^ I system the Jastrow method describes the wavefunction in terms of the product of two-body 
correlations between all pairs of nucleons, acting upon a suitable reference state. This ap- 
proach has been very fruitful in the description of extended strongly interacting systems, 
such as helium liquids at zero temperature and even nuclear and neutron matter. 

The practice of the Jastrow correlated method in infinite systems requires the use of 
massive resummation techniques, such as the hypernetted-chain (HNC) algorithms |^,^ or 
Monte Carlo (MC) sampling procedures . After more than 20 years of practice this method 
has proven to be perhaps the best variational description of many-particle systems. 

A different question is the treatment of finite strongly interacting systems, the main 
difficulty being related to the A-hodj interparticle correlation induced by the localization of 
the system. From the point of view of massive resummation of diagrams a la HNC, the most 
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unpleasant fact is that any diagram becomes elementary, which in turn makes the HNC 
resummation very involved. Even though the formal equations were written many years 
ago [0, it is only quite recently that some results for finite nuclei with model, realistic-like, 
interactions have been presented ^j. On the other hand, there is no special difficulty, not 
present in the case of extended systems, with regard to the use of variational MC methods. 
In any case, the presence of state-dependent interactions and, in turn, of state-dependent 
correlations gives rise to a formidable computational complexity. 

An alternative way to calculate expectation values for Jastrow wavefunctions consists 
of the consideration of cluster expansions, in which the mean value of an operator for an 
yl-particle system is obtained by means of nonlinear expansions of the expectation value for 
n-body subsystems. As in the case of HNC and MC methods, the basic idea is taken from 
statistical mechanics. Cluster expansions adapted to finite systems were analysed by Clark 
and Westhaus and used in nuclei by one of us [0. MC studies 0], as well as the recent 
HNC analysis have confirmed the quality of these expansions. 

Previous calculations of the binding energy for several light nuclei make use of state- 
independent (SI) Jastrow correlations between the nucleons. Our aim is to study the effect of 
dealing with state- dependent (SD) correlation functions within a Jastrow variational context. 
As an initial step we will concentrate in the simple case of ^He, as well as in simple central 
(but spin-isospin dependent) forces of semirealistic nature. Such interactions are usually 
termed as V4, the number 4 referring to the four two-body operators {U, <?« ■ a^, Ti ■ tj, {di ■ 
(Tj){Ti ■ Tj)}. The main difference with respect to variational MC studies is that we will carry 
out all calculations by means of algebraic techniques, avoiding the presence of statistical 
errors. This will be very important when dealing with heavier nuclei. 

The trial wavefunction, formally written as = FA\^), is the product of a Jastrow 
correlation factor F and an uncorrelated shell-model state. The latter is constructed with 
single-particle orbitals from the harmonic-oscillator (HO) potential, in order to be able 
to exactly remove the spurious centre-of-mass (CM) motion. Denoting by {(f)k{xi), k = 
1, . . . ,A} the set of relevant normalized single-particle states (specified for the coordinates 
Xi = {fi,ai,fi) of the ith nucleon), the antisymmetric reference state is nothing but the 
Slater determinant A\^) = det{(j)k{xi))/ VA^.. The Jastrow factor is the symmetrized PJlO[| 



product of two-body correlation functions, also of V4 structure, for all pairs of nucleons. The 
number of factors that build up the total F is therefore Nc = A{A — l)/2. The correlation 
corresponding to a given pair (ij) is expressed as 

fix,, X,) ^ h, = rinj) i + r (r,,) P^itj) + rin,) p,{tj) + nr,,) PMj)PAtj) (i) 

where the spin and isospin dependence is explicitly written in terms of spin or isospin 
exchange operators, and f^, , and /°^^ are scalar functions of the interparticle distance 
Tij = \ fi — fj\. The SD correlated g.s. wavefunction a la Jastrow is finally written, up to a 
missing normalization factor, as 

^(xi,...,a;^) = n /(xi,a;j)|^<l>(a;i,...,a;^) (2) 
where S is the symmetrizer for the operator inside the curly brackets. 
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It is worth stressing that the antisymmetrization of |$) — which is required to fulfil 
the fermionic nature of the nucleons — as well as the symmetrization in F — imposed by 
the (in general) non-commutative character of the different pair-correlation operators — are 
sources of notable technical complications. From the symmetrization of /12/13 • • • Ja-iAi we 
see that Nc\ terms indeed constitute the correct, complete correlation. Then, as long as 
the number of particles is small the problem remains manageable but its difficulty rapidly 
grows with increasing A. In passing, let us mention that in most works posing this sort of 
treatment |jlO|-Q , the so-called independent-pair approximation^ has widely been applied 



and so the question of symmetrization has been avoided. 

In order to obtain the g.s. energy we should compute the mean value of the Hamiltonian 

^ = -|^EV.^+ E V{x.,x,) (3) 

i=l j>i=l 

where /m = 41.47 MeV fm^, and the potential V is a two-body operator that depends 
on the relative spatial coordinate and the spin-isospin degrees of freedom of each pair of 
nucleons. The interactions employed here have V4 structure, which will be transformed to 
an exchange-operator basis as in the case of the correlation (|I]), 
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V{xi, Xj) = J2 9u{rij) J2 ^«.p Op{ij) (4) 

where Oi(^j) = i, Oalu) = Pa{i3). 0^]) = Priij) and 6^]) = Pa{i3)Pr{ij). This 
representation is very economical because in matrix form any exchange corresponds to a 



sparse matrix with a single non-null element per row and column |13[. For our calculations 
we will consider the quasirealistic potentials proposed by Afnan and Tang [0,0 and by 
Malfliet and Tjon |T^, whose various components have Gauss and Yukawa spatial depen- 
dence respectively, with strengths Vu,p and corresponding form factors Quij") = exp(— /i„r^) 
and Quif) = exp(— /Xtir)/r. The Coulomb interaction will not be taken into account in any 
case. Notice that the Hamiltonian (Q) is symmetric under particle permutation and so is 
the operator F^HF . This property enables us to find the total kinetic energy as A times 
the corresponding value for a single particle, and the total potential energy as A^^^; times the 
contribution of a single pair of nucleons. 

The ^He results presented here can be obtained exactly by a direct evaluation and further 
minimization of the mean value of the Hamiltonian with a trial state of the form (^. 
Nevertheless, we have followed a cluster representation with the purpose of gaining some 
insight into this approach for future application to more complicated situations, namely 
for heavier nuclei or more realistic interactions. In that sense, the present example serves 
as a test for the convergence of the cluster expansion when SD correlations are taken into 



^The independent-pair approximation involves rewriting each fij as a scalar function of r^j, times 
an SD operator of the form U + Uij, and then retaining from the expansion of ni<j(l + ^ij) only 
those terms with labels that correspond to independent pairs of nucleons. 
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account. Besides, the approximation procedure is useful for seeking the region where the 
minimizing parameters of the correlation functions should be found. The methodology we 
followed to compute the required expectation values is the one developed in that we 
summarize below. 

Consider an n-particle subsystem {1 < n < A) whose Hamiltonian operator if„ has the 
form d^) with the replacement A ^ n. The n-body integrals are introduced as 

{n I I ?i I n! 

n Uf e^^" 5 n 4 E ep ■ ■ ■ ^kA^pn)) (5) 
i'<j' J [i<j J p=i 

where the first sums extend over all the single-particle states relevant for the given nu- 
cleus, and the last one runs over all permutations of particle indices (ep are their asso- 
ciated parities). The matrix elements involve, besides spatial integrations over the whole 
space, summations for the spin and isospin variables of the n nucleons. If one deals with 
spin/isospin saturated nuclei, these summations translate into corresponding traces of prod- 
ucts of exchange operators. Their evaluation, although quite long, can be made once at 
the very beginning of the calculations and the resulting values can be stored for further 
use (a detailed explanation on this subject can be found in |jl3|). In the above equation, 
symmetrization of an n-hodj correlation is required and therefore nc\ = {n{n — l)/2)! re- 
orderings are obtained (both on the left and on the right), in contrast with the NqI terms 
arising from the complete A-particle problem. The quantity we are interested in, the min- 
imum value acquired by (\I'|i/|^)/(\E'|\E') — Eqm, can be derived from the A-hodj integral 
since 



A=0 



Ja{0) ^ ' 



The second term in the g.s. energy takes proper account of the CM movement and is given 
by the kinetic energy of a particle in the lowest HO state, i.e. E^m = 3 /i^Q!^/(4m), a being 
the HO constant. 

The binding energy should be obtained by scanning the space of trial correlation func- 
tions. Except for light nuclei, this is impracticable with present-day computers because of 
the huge number of evaluations to be performed. It proves useful to introduce a multiplica- 
tive approximation procedure by means of the cluster integrals Jy which are defined by the 
recursion 

Jn = n ^-^^^ (7) 

for all n = 1, . . . , A. The interpretation of this chain of identities is quite straightforward: 
one can obtain the n-body integral as a product of different contributions arising from fewer- 
body interacting subsystems. This idea allows us to decompose the energy expectation value 
as 



( 




(8) 

A=0 
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For a nuclear system, one may assume that a reduced number, say n, of cluster integrals 
would suffice in determining J a with great accuracy, the remaining cluster integrals being 
close to unity so as to produce negligible contributions SE^'^-'"'^ to the energy (||). The 
A-hodj integral calculated up to the nth order turns out to be 

Ja^Ji^... jP = Jf ^ (9) 

with J7i,...,j7„ deduced from the recursive formulae (|^. Starting with n = 1, one can 
systematically construct a series of approximations for the energy by taking the logarith- 
mic derivatives of J^"\ A detailed analysis of the characteristics and convergence of this 
multiplicative cluster expansion employing SI correlations has been carried out in previous 
works Up to ffist order the energy has purely kinetic origin and amounts to 



2^,2 



h a 



m 



3 + 5 



A 

T 



(10) 



This expression is valid for s- and p-shell nuclei, namely ^He, ^Be, ^^C and ^^O, for which 
the spatial single-particle states are 0s(r) = (a/ a/tt)^''^ exp(— and {(p^, (j)y, <Pz}p{^ = 
^/2a {x, y, z}(l)s{r). The ffist non-trivial result is obtained with one- and two-body clusters: 



— — 2 

J2 Jl 



(11) 

A=0 



Higher orders are constructed following similar steps, i.e. calculating the successive correc- 
tions as suggested in equation (||) so that 

^(n) _ ^(n-l) ^ 1^^^ ^^(n)_ (^3) 

Let us mention that at each order one can add at once the kinetic contribution of only one 
of the n particles and the potential due to a given pair of them, with corresponding factors 
n and ric- The approximation procedure will be performed up to a maximum of a few 
clusters, studying the convergence of the partial results. Note that very precise calculations 
are required, for two main reasons: the cancellations in 6E'^'^'> implied by the definition of 
Jn (as the n-body integral divided by all the Ju<n)i and the fact that the statistical factors 
(exphcitly factored out in the preceding equation) introduce an enormous scaling with n. 

Let us now specify the form of the pair correlation functions. A very useful, and rather 
easy to handle, parametrization is an expansion in a set of Gaussians. Our ansatz is the SD 
linear combination 

4- = i + E e"^'"^- E Op{tj). (13) 

m=l p=l 

The correlation lengths Pm and correlation depths am,p, together with the HO constant 
a, constitute a set of (real) free parameters to be determined from energy minimization 
(throughout this work, distances are measured in fm and the variables (3^ and a will be 



5 



given in fm~^ and fm~^ respectively). As it is clearly seen, the SI problem arises as a 
particular case by setting Om^p^i = and freely determining the 2Np + 1 variables 0^,1 
and a. 

In the case of ^He nuclei, the problem substantially simplifies due to the fact that the 
reference state exhibits spatial symmetry. Indeed, the uncorrelated g.s. is made up of 
the four particles in the same spatial single-particle state, (f)s- As is known, the opera- 
tor Pf{ij)Pfj{ij)PT-{ij) characterizes the exchange of all coordinates of the nucleons i and 
j. Since Pr{ij) leaves the helium wavef unction unaltered, the spin and isospin exchange 
operators acting on ^<l>(xi, 0:2, X3, X4) give a change of sign, and consequently the action of 
Po- is equivalent to minus that of P^. The same holds in second and third order of the cluster 
expansion for ^He. Then, for this nucleus we are able to rewrite the correlation factors with 
a simpler state dependence, the ansatz (|1^) becoming 



4- = i + E e"^™'^^ [am,c i + am,.PMj)] (14) 



with central scalar and spin-exchange constituents only. The state-independent study is 
performed, as usual, by setting 0^ 0- = for all m. 

In tables | and ^ we present the results corresponding to the g.s. energy of the ^He nu- 
cleus, obtained using different trial wavefunctions and nucleon-nucleon interactions. Table | 
corresponds to the Afnan-Tang S3 potential [|T^ and its modified version MS3 while 
in table || the Malfliet-Tjon interactions MT I/III and MT V are considered. The 
four-particle system has been characterized by a Jastrow prescription with state-dependent 
correlations containing one and two Gaussian components. For the sake of comparison, we 
also show in the tables the concomitant state-independent energies (the SI values have ear- 
lier been computed in assuming /m = 41.50 MeV fm^). Together with the minima 



obtained in each case, we exhibit the optimum values for the HO constant a, correlation 
lengths (3m and correlation depths 0^,0 and 0^,0-, for m < A^^ with Np = 1 and 2. 

The S3 potential is parametrized in terms of five ranges fiu and has only singlet-even 
and triplet-even channels; the MS3 interaction incorporates the odd ones. Although for 
helium both potentials yield the same results employing SI pair correlation factors, the 
introduction of a dependence upon the spin of the nucleons distinguishes between those 
interactions. The MS3 shift with respect to the SI binding energy turns out to be somewhat 
smaller in magnitude than the S3 one: 3.7% (9.1%) against 3.9% (10.0%) respectively, for 
one (two) correlation lengths. This difference is related to the saturation properties of the 
modified Afnan-Tang interaction. 

The potentials proposed by Malfliet and Tjon are a superposition of an attractive and 
a repulsive Yukawa component. In the MT I/III case, the spin-singlet and spin-triplet even 
channels are split; for the MT V case it is assumed that these two forces can be replaced by 
an average effective potential which is identical in both channels. The last situation allows 
one to treat the four nucleons as identical spinless bosons, meaning that the SI and SD 
approaches for the MT V potential give the same result, which is also consistent with the 
state-independent treatment for the MT I/III force. We find that with this interaction the 
state-dependent binding energy diminishes by almost 1% (3.7%) of the SI value for A"^ = 1 
(Np = 2). 
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The results presented here can be compared with other existing calculations, among 
which we may mention the energies given using MC techniques and the spin-isospin- 
dependent translationally invariant coupled cluster (TICI2) treatment |T^. In the case 



of the MT V potential, for instance, our [Np = 2)-result is very close to the Green function 
MC value of (-31.3 ± 0.2) MeV |g and the diffusion MC energy, (-31.32 ± 0.02) MeV 
T7| . After performing a coupled cluster calculation for the alpha particle, Zabolitzky |19 



arrived at the conclusion that the g.s. eigenvalue of the MT I/III Hamiltonian should be 
(—33.4 ± 0.1) MeV. The TICI2 method with SD linear correlations of the V4 type yields 
the following results fT^ for the quoted potentials: -28.19 MeV (S3), -27.99 MeV (MS3), 
-30.81 MeV (MT I/III) and -29.45 MeV (MT V). These values have not been optimized 
with respect to the HO parameter and hence they may be slightly underestimated in magni- 
tude. It can be seen that the results given here are comparable with the binding energies of 
the TICI2 method; moreover, our results with two betas are lower than the latter by more 
than 1 MeV. 

In order to look at the convergence of the above mentioned multiplicative cluster ex- 
pansion, we show the partial results corresponding to the minimum energies found with one 
and two (3 parameters. Table |T| contains for N/3 = 1 the first-order energy ^(1) as well 
as the concomitant corrections up to the second, third and fourth orders, which should be 
multiplied by the statistical factors 6, 4 and 1 respectively. One can see that in this situa- 
tion the SI fourth-order corrections are much smaller than the respective third-order terms, 
and using SD correlation functions |(5£'*^^^| remains lower than |(5£'''^^|; finally, the third- and 
fourth-order corrections have opposite signs. The influence of the introduction of a second 
SD correlation component is shown in table |V|, giving evidence that in this case four-body 
SD contributions gain in importance and thus are not negligible. Moreover, for the MT 
potentials the third- and fourth-order corrections both contribute to lower the second-order 
energy. We stress that, at this step, our purpose is to compare the minimum total en- 
ergy with the third-order approximation obtained assuming the optimum parametrization 
of tables | and ||, which in principle should not be equal to the minimum 

In conclusion, let us remark on the good agreement of our results with the available 
MC ones, when we take into account just two correlation components. We have shown the 
importance of introducing a dependence with the spin and isospin of the nucleons in the 
pair correlation functions, the shift in the energies with respect to the state-independent 
values being of around 1-10%, depending on the two-body interaction considered. As a 
side benefit, in this simple case of helium nuclei we were able to test the convergence of the 
multiplicative cluster expansion, a crucial treatment to be taken into account for the study 
of heavier systems because the 'size' of the problem is considerably reduced with respect 
to the full energy computation. We note that, at least in the case of ^He, the fourth-order 
contribution to the binding energy is non-negligible for the (SD, Nf^ = 2) approach. 



This work was partially supported by DGICyT (Spain) grant PB92-0820. M.P. acknowl- 
edges the Instituto de Cooperacion Iberoamericana (ICI, Spain) for a fellowship. 



7 



REFERENCES 



[1] S. Fantoni and S. Rosati, Nuovo Cimento A 25 (1975) 593 
[2] E. Krotscheck, Nucl. Phys. A 317 (1979) 149 
[3] W.L. McMillan, Phys. Rev. A 138 (1965) 442 

[4] S. Fantoni and S. Rosati, Nucl. Phys. A 328 (1979) 478 

[5] G. Co', A. Fabrocini, S. Fantoni and I.E. Lagaris, Nucl. Phys. A 549 (1992) 439 

F. Arias de Saavedra, G. Co', A. Fabrocini and S. Fantoni, Nucl. Phys. A 605 (1996) 359 

[6] J. Clark and P. Westhaus, J. Math. Phys. 9 (1968) 131 
P. Westhaus and J. Clark, J. Math. Phys. 9 (1968) 149 

[7] R. Guardiola, Nucl. Phys. A 328 (1979) 490 

[8] M.C. Bosca, E. Bucndfa and R. Guardiola, Phys. Lett. 198B (1987) 312 
[9] V.R. Pandharipande and R.B. Wiringa, Nucl. Phys. A 266 (1976) 269 

[10] J. Lomnitz-Adler and V.R. Pandharipande, Nucl. Phys. A 342 (1980) 404 

[11] J.C. Owen, Ann. Phys. NY. 118 (1979) 373 

[12] J. Carlson and M.H. Kalos, Phys. Rev. C 32 (1985) 2105 

[13] R. Guardiola, P.I. Moliner, J. Navarro, R.F. Bishop, A. Puente and Niels R. Walet, Nucl. Phys. 

A 609 (1996) 218 
[14] I.R. Afnan and Y.C. Tang, Phys. Rev. 175 (1968) 1337 

[15] R. Guardiola, A. Faessler, H. Miither and A. Polls, Nucl. Phys. A 371 (1981) 79 
[16] R.A. Malfliet and J.A. Tjon, Nucl. Phys. A 127 (1969) 161 

[17] R.F. Bishop, E. Buendfa, M.F. Flynn and R. Guardiola, J. Phys. G: Nucl. Part. Phys. 18 
(1992) L21 

[18] J.G. Zabolitzky and M.H. Kalos, Nucl. Phys. A 356 (1981) 114 
[19] J.G. Zabolitzky, Phys. Lett. lOOB (1981) 5 



8 



TABLES 



Potential 


Correlation 


Energy 


a 








S3, MS3 


SI, Nfi = 1 


-24.4042 


0.8293 


2.0502 


-0.7191 


- 


S3 


SD, Np = 1 


-25.3598 


0.8380 


1.9433 


-0.7012 


0.1128 


MS3 


SD, = 1 


-25.3119 


0.8373 


1.9486 


-0.7023 


0.1073 


S3, MS3 


SI, Np = 2 


-27.2136 


0.5795 


1) 1.6206 


-1.5596 












2) 0.2402 


1.0872 




S3 


SD, = 2 


-29.9378 


0.6617 


1) 1.4646 


-1.4704 


-0.2805 










2) 0.3892 


0.9511 


0.4877 


MS3 


SD, Nf) = 2 


-29.7034 


0.6570 


1) 1.4703 


-1.4617 


-0.2454 










2) 0.3765 


0.9413 


0.4409 



TABLE I. For the Afnan-Tang potentials, g.s. energy of a ^He nucleus and optimum correlation 
parameters for different situations (energies are given in MeV, a in fm~^ and f5m in fm~^). 



Potential 


Correlation 


Energy 


a 


f3m 






MT I /III, MT V 


SI, = 1 


-29.0604 


0.8174 


3.8314 


-0.7934 




MT I /III 


SD, Ns = 1 


-29.3460 


0.8201 


3.7559 


-0.7853 


-0.0642 


MT I /III, MT V 


SI, = 2 


-30.8752 


0.5089 


1) 3.4460 


-1.9090 












2) 0.1455 


1.4125 




MT I /III 


SD, Np = 2 


-32.0107 


0.6077 


1) 3.3441 


-1.4139 


0.2230 










2) 0.2018 


0.8203 


-0.3998 



TABLE II. Same considerations of table I, for the potentials proposed by Malfliet and Tjon. 
For the MT V interaction, being of Wigner type, the SD approach gives the same binding energy 
as the SI one, i.e. the optimum am,a are zero. 



Potential 


Correlation 




5E(^^ 


5E(^^ 




S3, MS3 


SI 


85.5600 


-15.5220 


1.1913 


-0.2076 


S3 


SD 


87.3595 


-16.0471 


1.4832 


-0.5300 


MS3 


SD 


87.2105 


-16.0198 


1.4785 


-0.5151 


MT I/III, MT V 


SI 


83.1312 


-15.6007 


0.5628 


-0.0559 


MT I/III 


SD 


83.6763 


-15.7429 


0.6180 


-0.1178 



TABLE III. Contributions (given in MeV), up to all cluster orders, to the minimum g.s. energy 
for both Afnan-Tang and Malfiiet-Tjon interactions and the optimal one-beta parametrizations. 
The CM corrections have not been included. 
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Potential 


Correlation 








,5^(4) 


S3, MS3 


SI 


41.7835 


-8.1422 


-2.4407 


0.0645 


S3 


SD 


54.4792 


-11.7120 


0.7360 


-3.4694 


MS3 


SD 


53.7015 


-11.5075 


0.5723 


-3.2238 


MT I/III, MT V 


SI 


32.2192 


-6.7370 


-3.6119 


-0.1697 


MT I/III 


SD 


45.9464 


-10.3470 


-0.0580 


-4.1564 



TABLE IV. Same considerations of table III, for the optimal two-beta parametrizations. 
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